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<* Abstract 

We study the spin excitations and the transverse susceptibility of 

a two-dimensional antiferromagnet doped with a small concentration of 

ON, 

Q\ | holes in the t-J model. The motion of holes generates a renormalization of 

-(— » ' 

d ' the magnetic properties. The Green's functions are calculated in the self- 

consistent Born approximation. It is shown that the long-wavelength spin 

G , waves are significantly softened and the shorter-wavelength spin waves 

O 

o 



become strongly damped as the doping increases. The spin wave veloc- 
ity is reduced by the coherent motion of holes, and not increased as has 



been claimed elsewhere. The transverse susceptibility is found to increase 
considerably with doping, also as a result of coherent hole motion. Our re- 
sults are in agreement with experimental data for the doped copper oxide 
superconductors. PACS numbers: 74.25.Ha, 75.40. Cx, 75.40Gb, 75.50.Ep 



1 Introduction 

After the discovery of high-T c superconductors there has been intensive inves- 
tigation of the magnetic properties 1 of doped copper oxide materials because 
of their connection to high temperature superconductivity. The undoped com- 
pounds are antiferromagnetic (AF) insulators. Doping introduces holes, 2,3 the 
charge carriers, in the AF square lattice of the copper oxide planes. The long- 
range AF order rapidly disappears at low doping, and superconductivity arises 
upon further doping. Strong two-dimensional AF fluctuations are nevertheless 
observed 4 even at fairly high doping, suggesting a conducting phase that, in 
spite of being paramagnetic, exhibits short-range AF order. A striking feature 
of the copper oxides is the strong sensitivity of their magnetic properties to 
the hole concentration, <5. Experiments have shown important softening and 
damping in the spin excitations, 5-7 as well as a significant increase in the spin 
susceptibility, 8 " 10 for the doped copper oxides. It is therefore important to 
study the interplay between doping and antiferromagnctism for an understand- 
ing of these materials. 

It is believed that the essential physics of strong electron correlations in the 
copper oxide planes is described by the t-J model 

<i,j>,cr <i,j> 

where c ia and Ci a are creation and annihilation electron operators acting on 
a reduced Hilbert space with no doubly occupied sites, the spin operator is 



Si = 5 Y,ap c \ a a ap c iPi n * = n *T + n a and n iv = c L c i"- In tnc C0 PP cr oxides, 
J ~ 1500^ and t ~ 3J. For the undoped materials, i.e. at half-filling, only 
the Heisenberg part of the Hamiltonian is relevant and it describes a spin- 1/2 
AF insulator. With doping, and nearly half-filling, the Hamiltonian describes 
holes moving in an AF background, the holes strongly interacting with the spin 
array. The motion of holes tends to disrupt the AF order, because a moving 
hole leaves behind a string of flipped spins. The charge carriers are then holes 
dressed by a cloud of spin excitations. 

The propagation of a single hole in a two dimensional antiferromagnet has 
been studied with a variety of approaches. Considering the t-J model in a 
Schwingcr boson representation, hole motion was treated within a self-consistent 
Born approximation (SCBA). 11-15 It was found that a hole can propagate co- 
herently because of its strong coupling to the spin excitations, having a quasi- 
particle bandwidth ~ J and energy minima at momenta q^ = (±7r/2, ±7r/2). 
The calculated spectral density shows a quasi-particle peak of intensity ~ (J/t) 2 ' 3 
and a broad incoherent multiple spin wave continuum, of width ~ 2zt (z is the 
coordination number), that is located at higher energies. These results are in 
good agreement with those from exact diagonalization of small clusters. 16 The 
study of hole motion has been extended for a finite concentration of holes, within 
the SCBA. 17 ^ 19 The results obtained show that, to first order in 5, the quasi- 
particle characteristics remain essentially the same, supporting a description of 
the quasi-holes as nonintcracting particles, filling up a Fermi surface consisting 
of four pockets located at momenta q^, having an enclosed area proportional 



to S. The spectral density for finite doping again contains a coherent and an 
incoherent part. The role of the coherent motion of holes on the magnetic prop- 
erties of the copper oxides has been a matter of discussion. Some authors, 17 ' 20 
in contradiction to others, 21 ~ 23 have claimed that the coherent motion of holes 
leads to stiffening of the spin excitations, while it is the incoherent motion of 
holes that generates significant softening, leading to an overall softening. 

In this work we study the effects of hole doping on the spin excitations, cal- 
culating both the softening and the damping, and determine the doping depen- 
dence of the transverse spin susceptibility of a two-dimensional antiferromagnet, 
discussing in particular the contributions from the coherent and the incoherent 
motion of holes. Our study starts from the t-J model in the Schwingcr boson 
representation and is carried out in the SCBA. It is shown that the magnetic 
properties are very sensitive to hole doping, as a result of the strong interaction 
between the hole and spin systems, and that the coherent motion of holes leads 
in fact to softening of the spin excitations. 

2 The Interaction Between Holes And Spin Waves 

Our system is described by the t-J Hamiltonian (1) on a two dimensional square 
lattice. In order to enforce no double occupancy of sites we use the slave-fermion 
Schwinger Boson representation Ci a = f}bi a , where the slave-fermion operator 
/• creates a hole and the boson operator bi a accounts for the spin, subject to 
the constraint fjfi + J2 a b la^ = 2S - 



We consider the low doping regime, S <C 1, where the states are close to 
the pure AF state, and hence exhibit long-range order. The AF state is ap- 
proximated by the Neel state, which in the Schwinger representation can be 
interpreted as a condensate of Bose fields b^ — V25 and bj[ = V2S, respec- 
tively in the up and down sub-lattices, and the bosons bii — bi and bj-\ = bj are 
then Holstcin-Primakov spin-wave operators on the Neel state. 

The Hamiltonian (1), with S = 1/2, then becomes 



Ht-j = -* E [h.f}{b\ + b 3 )+H.c. 



b\bi + b)b 3 + b % b 3 + b\b) 1 



(2) 



The transfer part describes the reversal of spins as the hole moves. In the 
Heisenberg part, the factor ( 1 — ///» ) ( 1 — flfj ) accounts for a loss of mag- 
netic energy due to doping, and for small hole concentrations it may be replaced 
by 1 - f}fi - f}fj. 

Applying Fourier transforms and the Bogoliubov-Valatin transformation for 
the spin variables: 6_ k = Uk/?l k + fkAc and 6k = t>k/3_k + u k/?k, where Uk = 
[((1 - 7kT 1/2 + 1) /2] 1/2 , vu = -sgn( 7k ) [((1 - 7k T 1/2 - l) /^'\ and 7k = 
\ (cosk x + cosky), we obtain from (2) the effective Hamiltonian 

H = -"/= E A/q-k [V(<h - k )/5-k + ^(q - k, k)/3 k ] + E^4/3k. (3) 

q,k k 

1 (1 

Here, V(q, k) = zt (7 q Uk + 7 q +kWk), w k = (zJ/2) (l — 7 k ) , the coordination 
number is z = 4, the sums run over the Brillouin zone for an antifcrromagnet on 



a square lattice, and N is the number of sites in each sub-lattice. In (3), the first 
term represents the interaction between holes and spin waves resulting from the 
motion of holes with emission and absorption of spin waves, and the second term 
describes spin waves in a pure antiferromagnet. In writing (3) we neglected an 
interaction term involving the scattering of holes by spin-waves, proportional to 
J, because its effect is small compared to the other term, proportional to t. 17 We 
note that at the bare level the holes have no dispersion. In fact, they propagate 
only after being dressed by spin waves. Here we study the renormalization of 
the magnetic properties induced by the dynamical interaction between the holes 
and the spin waves. 

3 Green's Functions For Spin Waves And Holes 

Given the magnitude of the couplings in the copper oxides, we make use of the 
Green's function formalism at zero-temperature. The Green's functions for the 
spin waves are defined as 

D-+(k,t-t') = -i(rp k (t)(3t(t')), 
D+-(k, t-t') = -i (r/£ k (t)/?_ k (f )) , 

D— (k, t-t') = -i (T0k(t)/3-k(f)) , 

D++(^t-t') = -i(T/3l k (t)dt(t')), 
where ( ) represents an average over the ground state. Their Fourier transforms 
satisfy the Dyson equations: 

d^ (k, lo) = d» v (k, W ) + J2 D ^ ( k > ^) n7 ' ( k > ") d6v ( k > w ) . 

7<5 



where /i, v = ±. The free Green's functions are 

D+- (k,«) = (-« - w£ + w?) _1 , 

D --(k,a;) = £>+ + (k, w ) = , 
with 77 — > + , and II 7 (k, u>) are the self-energies generated by the interaction 
between holes and spin waves. 

We calculate the self-energies in the SCBA, corresponding to only "bub- 
ble" diagrams with dressed hole propagators. These diagrams describe the 
decay of spin- waves into "particle-hole" pairs. The approximation neglects 
corrections to the hole-spin interaction vertex, which have been shown to be 
unimportant. 12 ' 15,17 The self-energies are given by 



1 /> + OC 7 

n 75 (k, W ) = -i-^C/^(k,q) / ^G(q,w q )G(q-k,w q -w) ) (4) 

where G(q, w q ) is the Fourier transform of the Green's function for the dressed 
holes, G(q,i - t') = -i (T/ q (t)/t(f )), and 



t/+-(k,q) = l/(q-k,k) 2 , t/-+(k,q) = V(q,-k) 2 , 
U— (k,q) = C/ ++ (k,q) = y(q,-k)y(q-k,k). 

The relations n^ + (k,w) = n +_ (-k, -w) and II — (k,w) = II ++ (k,w) are ver- 
ified. 

In the SCBA the holes are dressed by pure AF spin waves. This approach 
implies a spectral function for the holes that is composed of a coherent quasi- 

7 



particle peak and an incoherent continuum, with the quasi-holes filling up a 
Fermi surface that consists of pockets located at q^ = (±7r/2,±7r/2), as men- 
tioned above. We shall take for the hole spectral function the approximate 
form 

p(q,w) = [p coh (q,w)+p incoh (q,u;)] T ± (q)9(±c), (5) 

with, Fermi surface F~ (q) = J2i=i^(lF — \q — qi\), F + {q) = 1 — ^ _ (q), Fermi 
momentum qF = \r^5, and 

p coh {q,u) = a ^-£ q ), 
p incoh (q, u>) = M(\w\ - zj/2)6(2zt + z J/2 - |w|). 

Here the energies are measured with respect to the Fermi level, and the quasi- 
particlc dispersion can, near the minima at q^, be written as e q = e m in + 
(q — qi) 2 /2m, with an effective mass m ~ 1/J (neglecting band anisotropy). 
The quasi-particle residue is ao ~ (J/t) 2 ' 3 , and the remaining spectral density 
appears in the incoherent continuum of width 2zt and height h ~ (1 — ao) /2zt, 
satisfying the sum rule J ckvp(q,ui) = 1. 

The spin wave self-energies (4) are obtained in terms of the hole spectral 
function (5) by 



ir*(k, W ) = 1 J2 Ul5 (^ q) [ y (<i> k ; w ) + r (q - k > k ; -«)] . ( 6 ) 

with 



TV 
q 



Y(q,-k;u;) = / dco' dw 

JO J -cc 



„p(qV)p(q-k,w") 



u> + lu" — lu' + ir\ 
From (5) and (6) follows that the self-energies will present three contributions, 

ir 5 (k, U ) = n£(k, u) + n;t(k, uj) + n7 c * jc (k, u), 

corresponding, respectively, to transitions of holes within the coherent band, 
between the coherent and incoherent bands, and within the incoherent band. 
We have calculated these different contributions to lowest order in the hole 
concentration 6. The imaginary parts of these contributions are non-zero only 
in certain regions of the (k, to) space: ImII c c ^ for [—kqp/m + k 2 /2ml < u> < 
[kq F /m + k 2 /2m] , ImII Ci j c ^ for z J/2 < lu < zJ/2 + 2zt, and ImIIj Cjlc ^ 
for zJ<lu<zJ + 4zt. 

4 Magnetic Properties 

We now present the calculation of the effects of hole doping on the magnetic 
properties. 

The renormalizcd spin wave energy w^ is given by the poles of the Green's 
functions D(k, w), determined by the condition 

[(AT+r 1 - n+-] [(ATr 1 - n-+] - n++n~ = o. 

In the region where ImII(k, lu) = 0, we find to lowest order in 5, 

w k = w£ + R*n + -(k,w£), (7) 



leading to 

c^ = w£(l-r(k)), (8) 



where 



t\ 2 (1 ( k 2 



r{k] = 6a l\-j) [i{T^) m ' lh -- k)+ 



9(k-2q F ) 



>m 2 ky\ (I -7k~ (fc/2) 






/-t (l-a ) 2 [, , ao , ^ , ,*M f 1 f A: 3 



ln2 + -^_ln[l + 4- 
1 - ao V J 



sin 2 k x + sin 2 k y 

In r(k) the first term is generated only by the coherent motion of holes, i.e. 
II+~, whereas the second involves the incoherent motion resulting from the 
sum nj"~ + n^,^ c . One finds that both the coherent and the incoherent motion 
of holes generate a reduction of the spin wave energy, and hence give rise to 
softening of the spin excitations. The fact that the coherent motion of holes 
leads to softening, even in the regime where the spin wave velocity is larger 
than the hole Fermi velocity is explained in detail in the Appendix. 
In the long- wavelength limit, k <C 1, one has 

Wk = ck, (9) 

with 

c = Z c c , (10) 
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where Co = zJ/(2V2) is the spin wave velocity for a pure antiferromagnet, and 
the renormalization factor is 

Z c = l-5al(j) . (11) 

For finite hole concentrations one has Z c < 1, which implies a reduction of the 
spin wave velocity with doping. This effect is generated only by the coherent 
motion of the holes, as can be seen from (8). 

In the region where the spin wave dispersion crosses the pair excitation 
continuum, defined as the region where ImII(k,w) =/= 0, one finds, to lowest order 
in 5, that the spin excitations become damped, acquiring an inverse lifetime 
given by 

T(k) = -2Imn+-(k,w k ) (12) 

One finds that the damping is determined only by the coherent motion of holes, 
i.e., Imn+~ , because the contributions involving the incoherent motion, Imllj~~ 
and ImII+~ c , vanish in the relevant region of the (k, u>) space. Hence we have 

r 2 



r(k) = zjVSat ^ j ^ k{l _ il)l/2 F+- (k) x (13) 

v/1 - S 2 (<?k)0(l - k(5k)|) - v/1 - s 2 H?k)0(l - |*(-0k)|)" , 



with 



F + (k) = (cos k y — cos k x ) {cos(g-^k x ) — cos(g-^k y )) + 

-2(1 - 7 £) 1/2 (sin k x sm{g^k x ) + sink y sm{g k k y )) + 4(1 - 7 £), 
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where s(g k ) = (1 - gu)k/2q F , g k = (2w k /J)/fc 2 , and w k is given by (8). We 
note the strong doping dependence of the damping ~ yS, as compared to that 
of the reduction of the spin wave velocity (Z c — 1) ~ S, 6 -C 1. 

From (13) one has that for sufficiently small doping, long-wavelength spin 
waves remain well defined, whereas the shorter wavelength spin waves are damped, 
decaying into " particle- hole" pairs. As the doping increases more spin waves, 
in the shorter wavelength side, dive into the pair excitation continuum and be- 
come damped. For hole concentrations above a certain threshold 5* , such that 
the spin wave velocity equals the Fermi velocity, Z*ca/(k F /m) — 1, the spin 
wave dispersion lies entirely in the pair excitation continuum, and then even 
the long-wavelength spin waves are damped. In the limit k <JC 1 one has r ~ k, 
which implies that the spin waves are overdamped. 

The transverse spin susceptibility is defined by 

X±=x±(k = 0,w = 0), (14) 

where the dynamical susceptibility is given by 

/>oo 

X±(k,«) =» / dte iujt < [S x (k,t),S x (-k,0)} > . 
Jo 

Writing the spin operator in terms of the electron creation and annihilation 

operators, Sf= (Sf + S^)/2 with S+ = cLc^ and S~ = cl\Ci^ using the 

Schwinger boson representation with the bose condensation associated with the 

Neel state, and performing the Bogoliubov-Valatin transformation, one finds 

that the susceptibility can be expressed in terms of the spin wave Green's func- 

12 



tions by 

X x =-limf- — J [RcL>+-(k,0)+RcD ++ (k,0)] . (15) 

In (15) we have approximated < f}fib\ >~ S < b\ >, and neglected a prefactor 
(1 — 5) 2 which is caused by dilution of the spin lattice by holes. To lowest order 
in S, the transverse susceptibility is given by 



lim — —. r- 

k^o z.J(l + 7k) 



1 " zJ(l - t) 1/2 ( Ren+ ^ k '°) + RcII++ ( k ' )) 



(16) 
One finds that only the coherent motion, i.e., n+jf , contributes to the suscepti- 
bility in (16), because the contributions involving the incoherent motion, 11^ 
and n^ , vanish in the limit k — > 0. We then obtain 

X ± =Z xX ° ± , (17) 

where x° = l/(2zJ) is the transverse susceptibility for a pure Heiscnberg anti- 
ferromagnet, and the renormalization factor is given by 



2 / I 



Z x = l + 45ai\jj . (18) 

From (18) one sees that the transverse susceptibility increases with hole doping, 
this effect being determined by the coherent motion of holes. 
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5 Results and Discussion 

We have considered a two-dimensional antiferromagnet doped with a small con- 
centration of mobile holes and calculated the renormalization of magnetic prop- 
erties induced by hole motion. In our calculation it is assumed that there is 
long-range AF order in the system. In real materials true long-range AF or- 
der disappears at rather low concentrations, e.g., S c ~ 0.02 for La2~sSrsCu04. 
However, experiments have revealed that above such concentrations, there are 
large AF correlated regions in the system corresponding to the size of the mag- 
netic correlation length £, scaling like £ ~ l/\/<5. 24 Those regions can in particu- 
lar sustain spin excitations with wavelengths up to the region size. One expects 
that the results that we derived when there is long-range order, still describe 
the physics on length scales less than £, with £ large, when long range order is 
broken. 

We find that the spin excitations are very sensitive to doping, with signif- 
icant softening and damping occurring as a result of hole motion. In the low 
momenta region, the reduction of the spin wave energy is mainly determined 
by the coherent motion of holes, while the contribution from the incoherent 
motion becomes more significant with increasing momenta. The spin wave ve- 
locity decreases due to the coherent motion of holes, which for t/J = 3 and a 
concentration <5 = 0.02 produces a renormalization factor Z c = 0.96, while a 
concentration S — 0.05 leads to Z c = 0.90. However, for momenta k around 
qF, the slope of the spin dispersion shows a much higher reduction, by a fac- 
tor 0.91 for a concentration 6 ~ 0.02, and a factor 0.78 for a concentration 

14 



S ~ 0.05, as a result of the coherent plus the incoherent motion of holes. For 
S ~ 0.02 there is little damping since only spin waves near the upper end of 
the spin wave spectrum lie inside the pair excitation continuum. For 5 ~ 0.05 
the spin dispersion dives partialy into the pair excitation continuum, so that 
excitations with k > Iqp are strongly damped. For concentrations above the 
threshold 5* ~ 0.17, where the spin wave velocity equals the Fermi velocity all 
the spin waves lie in the pair excitation continuum, and therefore are completely 
damped. This occurs at a concentration well below the value for which the spin 
wave velocity would vanish. One expects that long-range order will collapse at 
a concentration S c < S* , implying a small value for the critical concentration, in 
agreement with experimental data. The disappearence of the long-range mag- 
netic order with doping will be discussed elsewhere 25 . Aeppli et al. 5 and Haydcn 
et al. investigated the spin dynamics of pure and doped LazsSrsCuOi, with 
5 = 0.05, and found that spin excitations within the AF correlated regions in the 
doped material show softening and damping with respect to the corresponding 
excitations in the pure material. Aeppli et al. found that the spin wave velocity 
in the doped material is renormalizcd by a factor 0.74(±0.08), while Hayden et 
al. found a renormalization factor 0.60. These values are to be compared with 
the renormalization factor for the slope of the spin dispersion around qF, the 
momenta range associated to the AF correlated regions, therefore our result, 
0.78 for 5 = 0.05, is in good agreement with experimental data. Experiments 7 
at a much higher concentration, 5 = 0.14, have also revealed a large broadening 
of the spin excitations with doping. 
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For the transverse spin susceptibility we find a significant increase with dop- 
ing, having for t/J — 3 a renormalization factor Z x = 1.17 for S = 0.02, and 
Z x = 1.42 for S = 0.05. This effect is also due to the coherent motion of holes. 
When long-range order is broken and the magnetic correlation length diverges, 
the susceptibility of the system should be essentialy given by x±. ■ An increase in 
the spin susceptibility with doping has in fact been observed experimentally, 8-10 
in agreement with our results. Above the critical doping a gap gradually opens 
in the spin excitation spectrum, and one may expect the magnetic correlation 
length in that regime to be determined by the imaginary part of the spin dis- 
persion. According to our results this implies an inverse correlation length 
proportional to Imll+J and therefore £ ~ 1/vS, precisely the scaling found 
experimentally. 24 

Spin excitations of a weakly doped antifcrromagnet have been investigated 
elsewhere in the SCBA. In Rcf.2I, the present authors considered the effect 
of the coherent motion of holes only, and showed that it generates spin wave 
softening, in agreement with the results in the present work. Other authors 17 ' 20 
claimed, however, that the coherent motion of holes leads instead to stiffening 
of spin waves. This contradicts our results, and may arise from approximations 
made in Refs.17 and 20. The renormalization of the spin excitations has also 
been calculated by another group, 23 but their calculation contains a self-energy 
independent of 5, a result difficult to understand. Becker and Mushelknautz in 
Ref. 22 also studied the effects of hole doping on the spin excitations, but using 
a different technique. They found softening and damping of the spin excitations, 
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due to both the coherent and the incoherent motion of holes, in agreement with 
our results. 

In conclusion, we have shown that the magnetic properties of a two-dimensional 
antiferromagnet are very sensitive to doping due to the strong interaction be- 
tween holes and spin waves, the coherent motion of holes leading to softening 
of the spin excitations, like the incoherent motion. 

Appendix 

The contribution of the coherent motion of holes for the rcnormalization of 
the spin waves excitations is given by Ren+~(k, w£). From (5) and (6) one has 

ReII+ c -(k,u;£) = al -^ 5>(|q-k|-(z F )% F -|q|)x 



><E 

where, e q = q 2 /2m. Given that 



C/+-(k,k-q- qi ) £/+-(k,q + q 4 ) 



Wk - ( £ q-k - £q) ^k + ( £ q-k - £q) 



^t/+-(k,k-q + qO=^C/ + -(k,q + qj )-(zt) 2 ^(7^ k+q ,-7^ +qi ), 



one has 



Rcn+- (k, w°) = a 2 ^ J2 (ll - k l - lF)e(q F - |q|) [A(k, q) - B(k, q)] , 



where, 



and 



A(k,q) = X>+-(k,q + qO T T^ ~ 



j=i 



K) - ( £ q-k - eq)' 
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4 l 

B(k,q) = (^) 2 ^(7q-k+ qi -7q+ q J r n , 



2 = 1 



K - ( £ q-k ~ £q)] ' 



For sufficiently small doping the spin wave velocity is larger than the Fermi 
velocity, and therefore A(k, q)0(|q — k| — qp)9(qF — |q|) > 0. However, one 
has that B(k, q) > A(k, q), even in the long-wavelength limit, k< 1, where 
B(k,q) - 2A(k,q). This implies that ReII+-(k,w£) < 0, and therefore soft- 
ening of the spin excitations, as given in (7). The inclusion of band anisotropy 
does not qualitatively change our results. 
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